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We develop a new method for calculation of quasi-normal modes of black holes, when the effective 
potential, which governs black hole perturbations, is known only numerically in some region near 
the black hole. This method can be applied to perturbations of a wide class of numerical black 
hole solutions. We apply it to the black holes in the Einstein-Aether theory, a theory where general 
relativity is coupled to a unit time-like vector field, in order to observe local Lorentz symmetry 
violation. We found that in the non-reduced Einstein-Aether theory, real oscillation frequency 
and damping rate of quasi-normal modes are larger than those of Schwarzschild black holes in the 
Einstein theory. 
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I. INTRODUCTION 

General Relativity is based on local Lorentz invariance, 
yet recently there appeared a lot of attempts to go be- 
yond local Lorentz symmetry Aether can be con- 
sidered as locally preferred state of rest at each point 
of space-time due-to some unknown physics. Einstein- 
Aether theory is general relativity coupled to a dynamical 
time- like vector field u a , which is called "aether". This 
theory is what comes instead of usual General Relativ- 
ity when local Lorentz symmetry is broken. Namely, u a 
breaks local boost invariance, while rotational symme- 
try in a preferred frame is preserved (see Q for a recent 
review). Different observable consequences of Einstein- 
Aether theory have been discussed recently [H . 

Quasi-normal modes have been studies recent years ex- 
tensively, because of their interpretation in Conformal 
Field Theory Q and possible, but yet unclear signifi- 
cance in Loop Quantum Gravity. At the same time, one 
of the most promising ways to check the theory of gravity 
as a fully non-linear theory is to observe the characteris- 
tic frequencies of black holes called quasi-normal modes 
QN frequencies are expected to be observed in the 
nearest future with the help of a new generation of grav- 
itational antennas. This suggests a unique opportunity 
to test the spontaneous breaking of the local Lorentz in- 
variance, through observing deviations of quasi-normal 
modes of black holes from their values in Einstein the- 
ory. 

Therefore, we would like to know what comes instead 
of the well-known quasi-normal spectrum of black holes 
in Einstein- Aether theory. For this one must have a solu- 
tion describing black holes in the Einstein- Aether theory. 
Fortunately such a solution has been recently obtained 
in Q, yet, only numerically. That was one of the mo- 
tivations for us for developing of the method for finding 
quasi-normal modes for potentials which are not known 



in analytical form, but are given only numerically in some 
region near a black hole (see Sec. HI of this paper). 
In case of asymptotically anti-de Sitter space-times, the 
method for finding QN modes for metrics given by a set 
of differential equations was proposed in |7|. 

We found here that quasi-normal modes for spheri- 
cally symmetric black holes in the non-reduced Einstein- 
Aether theory have larger real oscillation frequency and 
damping rate than the corresponding QN modes in the 
locally Lorentz invariant Einstein gravity. 

II. BASIC FORMULAS 

The lagrangian of the full Einstein- Aether theory forms 
the most general diffeomorphism invariant action of the 
space-time metric g a b and the aether field u a involving 
no more than two derivatives given by 

L = -R- K ab mn V a u m V b u n - \{9abu a u b - 1), (1) 

here R is the Ricci scalar, A is a Lagrange multiplier 
which provides the unit time- like constraint, 

K ab mn = c l9 ab g mn + c 2 5^5 b n + c 3 S^S b n + Ci u a u b g mn , 

where Ci are dimensionless constants. 

Spherical symmetry allows to fix C4 = 0. In this letter, 
following [||, we shall consider the so-called non-reduced 
Einstein- Aether theory, for which C3 = 0, and we can use 
the field redefinition that fixes the coefficient C2 [a]: 

c i 

C2 = ~2-4 Cl +3c?' 

so that ci is the free parameter. 

The metric for a spherically symmetric static black 
hole in Eddington-Finkelstein coordinates can be writ- 
ten in the form Q: 

ds 2 = N{r)dv 2 - 2B(r)dvdr - r 2 cft) 2 , (2) 



'Electronic address: 
t Electronic address: 



konoplya@fma. if . usp. br 
zhidcnko@fma. if. usp. br 



where the functions N(r) and B(r) are given by numeri- 
cal integration near the black hole event horizon One 
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can re-write this metric in a Schwarzschild like form: 

B 2 {r) 



ds z = -N{r)dr + 



N{r) 



dr z 



(3) 



Now, let us consider test scalar and electromagnetic 
fields in the background of the black hole given by the 
metric ([3]). Here we shall use general covariant gener- 
alizations of scalar and electromagnetic wave equations, 
i.e. we neglect interaction of these fields with aether. 
This can be well understood, since the background of a 
large astrophysical black hole has much more influence on 
propagation of fields during quasi-normal ringing than 
aether, at least not very far from the black hole. Yet, 
such phenomena as tales at asymptotically late times are 
probably affected by asymptotic regions of space-times 
as well as by a black hole itself. 

Thus, the wave equations for test scalar $ and electro- 
magnetic fields are: 

=0, (4) 



- A a ,,)g a VV=S)., = 0. (5) 

After making use of the metric coefficients ((3]), we found 
that the perturbation equations can be reduced to the 
wave-like form for the scalar and electromagnetic wave 
functions ty s and \& e /: 



d 2 * 



- + (w 2 - Vi(r))Vi = 0, 



B(r) 

drl ' v ~ ' H,// *' N(r) 
The effective potentials take the form: 

£(£ + !) 1 d N(r) 



V s = N{r)- 



r dr* B(r) ' 



dr. (6) 



(7) 



(8) 



Below we shall also use the particular case of the above 
metric when B(r) — 1 and N(r) = 1 — 2M/r, what cor- 
responds to an ordinary Schwarzschild black hole. Now 
we are in position to find quasi-normal spectrum for the 
above wave equations. Yet, as the function B(r), and 
N(r) are known only numerically, in the next section 
we shall have to develope a new technique for finding of 
quasi-normal modes in this case. 



III. A NEW METHOD FOR CALCULATION OF 
QUASI-NORMAL MODES FOR POTENTIALS 
UNKNOWN ANALYTICALLY 

For a wide class of problems of astrophysical interest, 
the dynamical wave equation has the form ^ , where the 
effective potential has the form of some potential barrier 
that approaches constant values at the event horizon and 



Figure 1: Potential for electromagnetic perturbations near 
Schwarzschild black hole (I — 2) and the same potential in- 
terpolated numerically near its maximum. Even despite the 
behavior of the two potentials are very different in the full re- 
gion of r, except for a small region near black hole, low-laying 
quasi-normal modes for both potentials are very close. 




spatial infinity. The general solution of the wave equation 
at infinity is 



* = A in ip in + A out ip ut, r* 



(9) 



The quasi-normal modes in this approach, by the defini- 
tion, are the poles of the reflection coefficients A out /Ai n . 
Out starting point is suggested by the WKB method Q 
where the asymptotic solutions of the wave equation near 
the event horizon and near spatial infinity are matched 
with Taylor expansion near the peak of the potential, i.e. 
between two turning points V(r) - co 2 =0. We state 
that the low laying quasi-normal modes are determined 
mainly by the behavior of the effective potential near its 
peak, while the behavior of the potential far from black 
hole is insignificant. Below we shall show that our state- 
ment is true. For this to happen, we shall consider the 
well-known potential for Schwarzschild black hole V(r) 
and also two other potentials which lay closely to the 
Schwarzschild potential near its maximum, but has very 
different behavior far from a black hole. These two po- 
tentials are chosen in the following way. We make a plot 
of the function V(r) of an analytic Schwarzschild poten- 
tial, then we find some number of points for this potential 
V ri near its maximum which serve us a basis for our first 
potential Vi n t which is an interpolation of these points 
near the maximum by cubic splines. Second potential 
Vfn is a fit of the above plot near the maximum by a 
ratio of polynomial functions. The possibility of making 
use of the numerical interpolation of the potential near 
its maximum gives us possibility of finding QN spectrum 
of fields near the solutions which are not known in ana- 
lytical form. 

Let us explain the interpolation technique more pre- 
cisely. Suppose we know some numerical solution for the 
metric. In other words, we are able to calculate the met- 
ric coefficients in a finite number of points with desired 
precision. Then, by using formulas (|7l8p . one can find 
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Table I: Schwarzschild fundamental quasi-normal modes of electromagnetic perturbations calculated by different methods: by 
using exact potential, using interpolation of the potential near its peak by 20 and 50 points per 2M unit, using fit of the metric 
functions by 10 and 100 points found with 1% precision. 



s=l 


Exact potential 


20 points interpol. 


50 points interpol. 


Fit (10 points, 1%) 


Fit (100 points, 1%) 


I — 1, WxvJdo 
£ = 1,WKB6 


U.4y 1 ( 4U-U. loDzlzl 

0.496383-0. 185274i 


U.4y looy-U. loOZl / 1 
0.512474-0. 174892i 


U.4yl ( oy-u.loozloi 
0.496625-0. 184945i 


U.4oyooU-U. Io4yyoi 
0.494266-0. 184069i 


n /ioiaq^ n i q^i qi i 
U.4ylOoO-U. loOloll 

0.496328-0. 185244i 


£ = 2,WKB3 
£ = 2,WKB6 


0.914262-0.190130i 
0.915187-0.190022i 


0.914225-0.190139i 
0.916133-0.189337i 


0.914261-0.190132i 
0.915134-0.190010i 


0.910125-0.188874i 
0.911038-0.188768i 


0.914155-0.190099i 
0.915080-0.189991i 


£ = 3,WKB3 
£ = 3.WKB6 


1.313467-0.191262i 
1.313797-0.191234i 


1.313441-0.19127H 
1.313923-0.191101i 


1.313466-0.191262i 
1.313799-0.19123H 


1.307435-0.189994i 
1.307761-0.189967i 


1. 313311-0. 191230i 
1. 313641-0. 191202i 


£ = 4,WKB3 
£ = 4,WKB6 


1.706036-0.191730i 
1.706190-0.191720i 


1.706016-0.191740i 
1.706205-0.191683i 


1.706036-0.191730i 
1.706190-0.191720i 


1.698158-0.190457i 
1.6983 10-0. 190448i 


1.705833-0.191698i 
1.705987-0.191688i 


£ = 5.WKB3 
£ = 5,WKB6 


2.095741-0.191968i 
2.095826-0.191963i 


2.095725-0.191978i 
2.095818-0.191953i 


2.095741-0.191968i 
2.095826-0.191963i 


2.086038-0. 190693i 
2.086121-0.190688i 


2.095492-0. 191936i 
2.095575-0.19193H 



Table II: The same as in the table U for scalar field. 



s=0 


Exact 


20 points int. 


50 points int. 


Fit (10 points,l%) 


Fit (100 points 1%) 


£ = 1,WKB3 
£ = 1,WKB6 


0.582228-0. 196003i 
0.585819-0.195523i 


0.582003-0.19580H 
0.589883-0.192359i 


0.582219-0.195995i 
0.585940-0. 195409i 


0.579382-0. 194684i 
0.582932-0.194211i 


0.582155-0.195970i 
0.585745-0. 195490i 


£ = 2,WKB3 
£ = 2,WKB6 


0.966422-0. 193610i 
0.967284-0. 193532i 


0.966365-0. 193592i 
0.968078-0. 192999i 


0.966420-0.193610i 
0.967284-0. 193532i 


0.961834-0.192317i 
0.962685-0. 192240i 


0.966304-0. 193577i 
0.967166-0.193500i 


£ = 3,WKB3 
£ = 3.WKB6 


1. 350412-0. 193024i 
1.350732-0.19300H 


1.350384-0.193029i 
1.350850-0.192884i 


1. 350411-0. 193025i 
1.350739-0.192995i 


1.344061-0.191738i 
1.344377-0.191715i 


1.350249-0. 192992i 
1.350569-0.192969i 


£ = 4,WKB3 
£ = 4,WKB6 


1.734680-0.192793i 
1. 734831-0. 192784i 


1. 734662-0. 192802i 
1. 734848-0. 192750i 


1.734679-0.192793i 
1. 734834-0. 192782i 


1.726554-0.191509i 
1.726704-0.191500i 


1. 734471-0. 192760i 
1. 734622-0. 192752i 


£ = 5,WKB3 
I = 5,WKB5 


2. 119141-0. 192678i 
2. 119224-0. 192674i 


2. 119127-0. 192688i 
2. 119219-0. 192665i 


2. 119140-0. 192678i 
2. 119224-0. 192673i 


2.109235-0.191395i 
2.109316-0.191391i 


2. 118886-0. 192646i 
2. 118969-0. 192641i 



the values of the potential at these points with the de- 
sired precision and any number of its derivatives. It is 
clear that the precision of the higher derivatives decrease 
very rapidly, so that, in order to find them with a good 
accuracy we have to find a lot of interpolation points 
very precisely. If the accuracy is not good enough, the 
calculations lead to some random values. Thus we feel 
that the potential has to be very smooth (and for phys- 
ically relevant problems it is) to use the WKB approxi- 
mation, because WKB method requires some number of 
the derivatives at the maximum of the effective potential. 
The 6-th order WKB formula reads 



i=6 



5>= 



1 

2' 



(10) 



where the correction terms of the i-th WKB order A, can 
be found in [§] and [nj, Q = V — uj 2 and Q l means 
the i-th derivative of Q at its maximum. For the WKB 
formula of sixth order one needs twelve derivatives of the 
potentials. 

Now let us see how the above programs of fit and 
interpolation work for the well- known case of the 
Schwarzschild black hole. In cases when we know the 
metric functions N(r) and B(r) numerically in some re- 
gion near the black hole using these functions in the for- 
mulas for potentials ([7][8]) gives us numerical values for 
the potential near the black hole. The accurate plot of 
the effective potential ([5]) £ = 2 is shown on Fig. IIIII 
together with interpolated potential Vi n t and fit poten- 



tial Vfu where the interpolation is made on the basis of 
75 points from r = 2M until r = 5M and fit is made 
on the basis of 100 points found with 1% precision. Fit 
functions are chosen as fractions of two polynomials: 



N(r) 



i=N N 
i=0 



-, B(r) = 



i=N B 

E 

i=0 



i=N E 



i=l 



which are substituted into equations ([7]) and ([5]). The 
numbers ./V/v and Ng determine the number of terms 
in the polynomials and are chosen in order to provide 
best convergence of the WKB series. Coefficients a\ , 

bj , a\ B \ b\ B ^ are determined by the fitting procedure. 
We see from the above plot, that at larger r the inter- 
polation potential Vi n t grows, while the precise function 
and fit function are decreasing. In the region near black 
hole r = (2M, 5M) the difference between the potentials 
is negligible. The interpolation was performed by cubic 
spline with the help of Mathematica®. In order to di- 
minish accumulating of an error, upon calculating of the 
first derivative of the interpolated potential, we had to 
interpolate also the obtained values for the first deriva- 
tive and all twelve consequent derivatives up to the sixth 
order. 

Let us look at the tables lllllTl The known 6th-order 
WKB quasi-normal frequencies are compared with those 
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obtained by using formula (fTOf for interpolation function 
Vi n t and for fit function Vfu- Within our approach we 
have two kinds of errors: error due-to WKB approxima- 
tion, and error due-to approximation of the potential by 
fit or interpolation. The error because of the WKB ap- 
proximation can be estimated simply by comparison of 
the QN values within different WKB orders. The initial 
deviations from the accurate potential due-to fit or in- 
terpolation induce larger errors for higher derivatives of 
the potential, and thereby, show themselves as a final er- 
ror in the computed QNMs. Therefore, if one gets good 
convergence of the WKB series, it means that the devi- 
ations of the approximated potential from the accurate 
one are small enough, so that induced errors in higher 
derivatives are small. Thus, the critcrium of goodness 
of a fit or interpolation is convergence of the WKB se- 
ries. An interpolation, to give a good WKB convergence, 
can be improved in the following way: one should in- 
crease the density of interpolation points, and, at the 
same time, increase the accuracy with which each point 
is given. Unlike interpolation, the fit approximation is 
more economic and works better with increasing either 
number of fit points or accuracy with which each point 
is given. 

Using the above mentioned procedure of improving of 
the fit and interpolation, we can, in principle, decrease 
the error due-to potential approximation until as small 
value as necessary. From tables [TilTTl we see that larger 
number of points for interpolation or fit gives better ac- 
curacy as compared with data for the exact potential. 
While for interpolation we had to increase precision of 
the potential values, for fitting it was enough to use the 
fixed precision of two digits after point for simple illus- 
trative Schwarzschild case. For instance, for t = 1, s = 1 
mode the difference between the QN frequency for the ex- 
act potential and for the fit potential is 0.011% for Reu> 
and 0.016% for Imu. The same order of error holds for 
higher multipoles. For lower multipoles, the error due-to 
fit approximation is two orders smaller than error due-to 
WKB technique. Therefore, when interpreting our new 
results in the next section, for Einstein- Aether theory, we 
may neglect error due-to potential approximation. 



IV. QUASI-NORMAL MODES FOR BLACK 
HOLES IN EINSTEIN- AETHER THEORY 

In the Einstein- Aether theory we know the metric func- 
tions N(r) and B(r) only numerically, as a result of some 
numerical integration in some region near the black hole 
Q. In order to integrate equation for metric coefficients 
with required precision, we had to find the accurate val- 
ues of /'(f/i) with the help of shooting method described 
in Q. Here we shall use the fit approximation, because, 
as explained above, it works better for this case. In this 
section we used for our fit 10000 points given with the 
six digit precision after the point. This is 100 times more 
points given with accuracy four orders better than that 



used in previous section for the Schwarzschild case. 

The polynomial fit for the metric or the effective poten- 
tial is characterized by the number of terms in its mul- 
tiplicator and denominator for metric coefficients N(r) 
and B(r), called here Njy and Nb respectively. There is 
some optimal number for which the convergence of WKB 
series is best. Practically we continue in the following 
way: in order to find optimal A/jv we search for minimal 
difference between third and sixth order WKB values, 
first for the case B(r) = 1. When we know an optimal 
value iVjv, in a similar fashion, i.e., by looking for best 
WKB convergence, we are in position to find the optimal 
fit for B(r). Quick WKB convergence shows that higher 
derivatives of the metric coefficients are calculated with 
the best accuracy. 

From tables IIIUlVl we can see that real and imaginary 
parts of uj is about 1% greater for c\ = 0.1. When c\ is 
growing up to c\ = 0.77, both Retu and Imu are grow- 
ing by about 10% and 20% respectively. For I > 1, the 
WKB method gives an error of about a fraction of a per- 
cent. Therefore our conclusion, that in the above region 
of ci < 0.77, the observed here effect of a few percents is 
much larger than the order of the WKB error. For the 
particular case of s = i = 0, the WKB error is about 5% 
(see Table [V}, therefore we cannot judge about the value 
of the aether effect in the case s — I = for small c\ . Yet 
for ci ~ 0.7, the effect a few times greater than the er- 
ror, and, qualitatively the behavior of the QN spectrum 
is the same: the larger C\ leads to a larger real oscillation 
frequency and damping rate. 



V. DISCUSSION 

The quasinormal modes of black holes should be fea- 
sible for observation by new generation of gravitational 
antennas. This would make a great impact on gravita- 
tional physics, where gravitational waves have not been 
not observed. 

At the same time, the gravitational consequences of 
Local Lorentz symmetry violation must show itself in ra- 
diative processes around black holes. Gravitational ra- 
diation damping of binary pulsars orbits reproduces the 
weak field gen eral relativity at the lowest post- Newto- 
nian order [TlJ. The significant difference between Ein- 
stein and Einstein- Aether theories can show itself in the 
regime of strong field, i.e. in changing the characteristic 
quasi-normal spectrum of black holes. Thereby, aether, if 
exists, could be indirectly observed through observation 
of characteristic spectrum of black holes. This motivated 
us to perform the present study of quasi-normal modes 
for black holes in Einstein- Aether theory. 

In this paper we developed a new method for finding 
quasi-normal modes for potentials determined only nu- 
merically in some region near a black hole. Then, the 
method has been applied here to numerical black hole 
solution in Einstein-Aether theory, which has been ob- 
tained recently in Q. 
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Table III: Fundamental QN frequencies of the non-reduced theory for s — 



Cl 


I = 1 


1 = 1 


£ = 3 


i = 4 


£ = 5 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.77 


0.589693 - 0.197691i 
0.593814 - 0.199635; 
0.598713 - 0.203470; 
0.604163 - 0.207262; 
0.610513 - 0.212271; 
0.618482 - 0.218866z 
0.628421 - 0.229250; 
0.637216 - 0.240653; 


0.973933 -0.195722; 
0.981227 -0.197457; 
0.989398 - 0.201414; 
0.998783 - 0.205259; 
1.009919 - 0.210209; 
1.023707 - 0.216876z 
1.042296 - 0.226654; 
1.060028 - 0.236792; 


1.360108-0.195176? 
1.370409 - 0.196886 i 
1.381943 - 0.200829; 
1.395215 - 0.204663; 
1.410990 - 0.209603; 
1.430563 - 0.216280; 
1.457064 - 0.226066; 
1.482451 - 0.236224; 


1.746918 - 0.194949; 
1.760201 - 0.196648i 
1.775082 - 0.200580i 
1.792215 - 0.204405; 
1.812598 - 0.209335; 
1.837918 - 0.216006; 
1.872255 - 0.225787; 
1.905212 - 0.235948; 


2.134014 - 0.194833; 
2.150271 - 0.196527; 
2.168489 - 0.200452; 
2.189472 - 0.204271; 
2.214446 - 0.209195; 
2.245488 - 0.215860; 
2.287621 - 0.225637; 
2.328101 - 0.235798; 



Table IV: Fundamental QN frequencies of the non-reduced theory for s = 1 



Cl 


I = 1 


1 = 1 


£ = 3 


I = 4 


£ = 5 


0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.77 


0.499500 - 0.187126z 
0.503184 - 0.188782; 
0.506385 - 0.192059; 
0.510456 - 0.195179; 
0.515081 - 0.199198; 
0.520683 - 0.204318; 
0.527472 - 0.211949; 
0.533464 - 0.219667; 


0.921357 -0.192128z 
0.928321 -0.193764; 
0.935562 -0.197618; 
0.944082 - 0.201334; 
0.954092 - 0.206114; 
0.966339 - 0.212544; 
0.982634 - 0.221933; 
0.997944 - 0.231610; 


1.322831 - 0.193366; 
1.332882 - 0.195021; 
1.343774 - 0.198918; 
1.356438 - 0.202690; 
1.371426 - 0.207550; 
1.389924 - 0.214116; 
1.414826 - 0.223723; 
1.438525 - 0.233670; 


1.718010 - 0.193859; 
1.731094 - 0.195524; 
1.745483 - 0.199429; 
1.762147 - 0.203217; 
1.781923 - 0.208800; 
1.806415 - 0.214705; 
1.839523 - 0.224382; 
1.871181 - 0.234419; 


2.110397 - 0.194105; 
2.126489 - 0.195775; 
2.144308 - 0.199683; 
2.164908 - 0.203477; 
2.189388 - 0.208370; 
2.219757 - 0.214992; 
2.260890 - 0.224700; 
2.300312 - 0.234780; 



We showed that quasi-normal modes of Einstein- 
Aether black holes have larger damping rate and real 
oscillation frequency than QNMs of Schwarzschild black 
holes have. We considered here scalar and electromag- 
netic test fields in the background of a spherically sym- 
metric black hole. The wave equation for gravitational 
perturbations could be found by consideration of pertur- 
bations of the complete system of Einstein- Aether equa- 
tions. Then our approach for finding of the quasi-normal 
modes can be applied as well. Yet, it is well known, that 
quasinormal modes in the eikonal approximation are the 
same for fields of different spin. Therefore we expect 
qualitatively the same results for gravitational perturba- 
tions as we observed for scalar and electromagnetic per- 
turbations. If the breaking of Lorentz symmetry is not 



very small, i.e. c± is not very small, the deviation of the 
QN modes from their Schwarzschild values might be ob- 
served in the forthcoming experiments with gravitational 
antennas. 
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